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Abstract 

We study the thermodynamic behavior of a decaying scalar field coupled to a 

relativistic simple fluid. It is shown that if the decay products are represented 

by a thermalized bath, its temperature evolution law requires naturally a new 

phenomenological coupling term. This "energy loss" term is the product 

between the enthalpy density of the thermalized bath and the decay width of 

the scalar field. We also argue that if the field (p decays "adiabatically" some 

thermodynamic properties of the fluid are preserved. In particular, for a field 

decaying into photons, the radiation entropy production rate is independent 

of the specific scalar field potential V{(p), and the energy density p and average 

number density of photons n scale as /3 ~ and n ^ . To illustrate these 

results, a new warm inflationary scenario with no slow roll is proposed. 
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In the new inflationary scenario accelerated expansion and reheating are separated into 
two distinguished periods. The flrst one is an exponential growth of the scale factor with the 
Universe evolving to a supercooled state. Due to this adiabatic expansion the temperature 
of the Universe decreases nearly 10^^ orders of magnitude At the end of this supercooling 
process the Universe is reheated. The fleld rapidly oscillates about the global minimum of 
its potential, and the energy density of the inflaton fleld is completely or almost completely 
converted into radiation in less than one expansion Hubble time. As a matter of fact, 
either on its early or modern version based on parametric resonance (sometimes called 
preheating) P,^, the reheating is a very fast and extremely nonadiabatic mechanism. 

In principle, if a sustained radiation component during inflation is allowed, the super- 
cooling and subsequent reheating could be supressed or at least weakened by many orders 
of magnitude |0. The flrst fleld motivated scenario based on this idea is the isothermal or 
warm inflationary picture as proposed by Berera [0. Like in new inflation, the warm picture 
starts from a high temperature phase transition with the universe evolving through a de 
Sitter inflationary period dominated by the scalar fleld potential. However, in the course 
of the expansion, energy is continuously drained from the fleld to the thermal bath. The 
whole process is described by an "energy loss" term, TcfP'. At the level of the scalar fleld 
equation of motion, it contributes like an additional viscosity r0, which may dominate the 
term "iHcj) corresponding to the redshift of (f) (the momentum of fleld) by the expansion. 
In the warm picture, the persistent thermal contact during inflation implies that the scalar 
fleld evolves in a sort of over damped regime. As a result, the Universe approaches a state 
where the dilution of the thermalized bath (due to expansion) may continuously be compen- 
sated by the production of particles from the decaying scalar fleld, thereby guaranteeing the 
constancy of the temperature. If this scenario works during an interval of time long enough, 
thermal fluctuations may also produce the primordial spectrum of density perturbations 
0J^(see also for an updated and more detailed analysis). 

More recently, this nice picture has severely been criticized in terms of its physical 
plausibility by Yokoyama and Linde 0]. The basic argument is that TcjP', is not sufficiently 



strong in the regime where it should describe warm inflation. This coupling seems to be only 
a small contribution in a sub-leading thermal correction, and as such, it is not expected to 
play a prominent role in the inflationary process |]lT|(see, however, Refs. and fT^]). On 
the other hand, as discussed in an extended framework [|12[, the basic idea of warm inflation 
is so attractive that it cannot be discarded without a more detailed investigation or even 
further attempts to solve the above pointed out difficulties. In particular, it is not so neat 
that an "energy loss" term like Tip"^ is a realistic approximation for describing the energy 
dissipated by the field to a thermalized bath. 

In what follows we apply basic thermodynamics arguments for studying the decay of the 
scalar field coupled to a thermalized bath which is represented by a relativistic simple fiuid. 
We recall that thermodynamics has often been used when the underlying microphysics of a 
given phenomenon has not been completely clarified. The basic reason is that it yields rela- 
tions among the macroscopic quantities whose validity is independent of the microphysics on 
which they ultimately depend. Therefore, it is interesting to consider some thermodynamic 
criteria on the question related to the plausibility of warm infiation. As we shall see, this 
macroscopic treatment clearly suggests a new form to the coupling term, which depends ex- 
plicitly on the created component through its enthalpy density. It has a quite simple form, 
namely, r0(p + p), where is the decay width of the scalar field, and may alter significantly 
the studies of the reheating period, as well as the original and extended warm infiationary 



pictures [13 



We will limit our analysis to homogeneous and isotropic FRW fiat universes. Following 
standard lines, the energy content is a mixture of a coupled scalar field plus a relativistic 
simple fiuid representing the thermalized bath. The total energy stress of this self-gravitating 
mixture obeys the Einstein field equations (EFE). In addition, since the scalar field works 
like a source of particles, the particle fiux of the relativistic bath (A^'^ = nu^) must satisfy 
a balance equation of the type N!^^ = \I'. The basic equations are: 



V{<P)+p] = 3H\ (1) 



5(y-^W+^) - -2^-^^- (2) 

n = —3nH + riT^, (3) 

where the dot means time derivative, H = ^ is the Hubble parameter, m^^ = 1/G is the 
Planck mass, p and p denotes the energy density and pressure, and n is the particle number 
density of the relativistic fluid which can be radiation though is not restricted just to that. 
The flrst quantity on the right-side represent the dilution of n due to the expansion of 
the universe, while the second describes phenomenologically the source of particles, that is, 
the rate produced by the fleld 0, which has been parameterized by the decay width of the 
scalar fleld (\[' = nF^). Note that can be rewritten in more enlightening form 

where = nR^ is the number of fluid particles per comoving volume. In particular, if 
is constant, the decaying of the scalar fleld increases N exponentially in the course 
of time. In the the lack of a detailed model, is supposed to be a generic function of 
the temperature. It should be noticed that and its possible dynamic influence has not 
previously been considered in the warm context. We argue here that such an approach is 
thermodynamically inconsistent because (|]) may change the evolution of the universe even 
in a single fluid description with no particle production [Q, and, indeed, it is fundamental 
for a consistent thermodynamic description of a decaying scalar fleld. Another point is 
related to the nature of the interaction term. The energy conservation law (m^T/"^.^), which 
is contained in the EFE, may also be written as a balance equation for the fluid energy 
density 

p + 3H{p + p) = -p^ -3H{p^ + p^), (5) 

or equivalently, 

p + 3Hip + p) = -u,T^\^. (6) 



and may substitute the second FRW equation. In warm inflation P,p^, as well as in the 



first studies of reheating |^ , the above equation was solved assuming the phenomenological 
coupling u^T^^ ^ = — r0^. However, although very useful due to its simplicity, a physically 
reasonable coupling term should be dependent on the nature of the created particles and not 
only of the scalar field itself. These two points are inextricably intertwined, in such a way 
that a proper macroscopic treatment (taking into account @) give rise to an alternative 
coupling term. 

The macroscopic quantities of the 7-fiuid are related by the Gibbs law 

nTda = dp — ^ ^ dn, (7) 
n 

where a is the specific entropy (per particle) and T is the temperature. 

By adopting T and n as basic thermodynamic variables, one may show from (^, (^) and 
(0) that the evolution temperature law is 



The first term on the RHS of (|]) is the usual equilibrium contribution. For an expanding 
fluid {R > 0), one finds T < 0, as it should be. The remaining terms display the out of 
equilibrium contributions due to the sources of particles (nF^) and energy {u^T^" ). A 
more convenient form to the temperature law is obtained inserting (^ into (|^) to obtain 
T f dp\ n 1 



T \dpj„n 

^ I dT 



(p + p)F0 + u^T, 



(9) 



The interaction term {u^Tj^'^ ^) can be determined only if one imposes some thermodynamic 
constraint in the above expression. In order to do a realistic physical choice, we first notice 
that the equilibrium thermodynamic relation for a radiation bath {p = |p), which is given 
by n oc T'^, is possible only if the second term on the RHS of (P) vanishes identically, that 



^mT;. = -(p + p)r<A- (10) 



In general, for a fluid satisfying the 7-law equation of state, p = (7 — l)p, from (|) the 
particle number scales as n oc T~ , as expected for a thermalized 7-bath, only if ([T0|) is 
satisfied. In addition, replacing (|1^) into (^, and using one may verify that p oc T~ . 
Since the whole argument may be inverted, it follows that the equilibrium relations for p 
and n are preserved, thereby characterizing a thermalized bath, if and only if the coupling 
term is proportional to the enthalpy density of the created particles as given by ([10|) . This 
coupling term is completely different from r0^, which has extensively been adopted in the 
literature. As should be expected in physical grounds, it depends explicitly on the nature 
of the created component (note the prefactor p + p = 7p on the RHS of (p!0|)). 

What about the underlying physical meaning of the new coupling term (p!0D ? Replacing 
(I) and (H) in the time derivative of (|^) one finds 

1 



riT 



(11) 



Therefore, if (|T0|) is valid, the specific entropy of the created particles remains constant. A 
similar condition has previously appeared in cosmologies with creation pressure [jl^], as well 
as in cosmologies with a phenomenological time-dependent A-term [|TB|, and was termed 
"adiabatic" creation, because the entropy and the total number of particles increases but 
the ratio o' = -j^ remains constant. The condition a = implies that the entropy 5* of 
the thermalized bath increases proportionally to the comoving number of created particles 
(A^ ~ nR^). More precisely, f = ^ = r</), which is a direct consequence of (|^). Naturally, 
even for "adiabatic" decaying, there are out of equilibrium contributions encoded in the 
temperature law (j^). Without loss of generality, let us analyze the case for a radiation bath. 
Inserting the coupling term ( p!OD into and (H), and rewriting (H) one finds the set of 
equations: 

^ - -- + (12) 

T ~ R ^ 3 ' ^ ^ 

4 

p^ + 3if(p0 + p^) = --Pr-r^, (13) 
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4 

pr + 4:Hpr = -pr^^ , (14) 

+ ^ = ^ ■ (15) 



SrirH 3H 

If the coupling term is negligible {T^ <C 3H), we see from ([T^ ) that the standard law is 
satisfied {TR = const), whereas (|T^) and (|T3p imply that Pr oc i?"^ and oc i?"'^ as 
should be for an adiabatic evolution. The opposite regime (T^ ^ 3H) defines an extreme 
theoretical situation, where the decay process is a phenomenon so fast that the dilution due 
to expansion is more than compensated (T > 0). An intermediary situation occurs if this 
ratio is of the order of unity [T^ ~ 3H). In this case, from ([I^ ) the bath temperature remains 
nearly constant (T ~ 0), and as a consistency check for the whole thermodynamic scheme, 
we see directly from (|1^) and (|15D that n and p are also nearly constant. Thus, a warm 



scenario with a constant radiation energy density is rigorously defined only when the scalar 
field decay rate is fine-tuned to the inverse of the Hubble time (F^ = 3H). This is exactly 
what one should expect from physical grounds: the dilution of the thermalized bath (due to 
expansion) is balanced by the particle production regardless of the dynamic details or any 
initial condition. Once this fine tuning is assumed, the dynamic behavior of the universe, if 
exponential (de Sitter), inflationary power law or even an ordinary decelerated FRW type 
expansion is determined by the energetics of the system. For example, if V{(j)) ^ <p'^ + pr, the 
resulting scenario resembles the exponential warm inflation originally suggested by Berera 
0. However, we see from (0) that the slow-roll conditions including the decay rate will be 
modified. We recall that the standard slow-rollover period requires F ^ if so that pr ^ 0^ 
(see Eqs.(13) and (17) of and also (7) of [|15|)- Within our approach these additional 



conditions are not necessary, because the thermodynamic 

constraint defining warm inflation is now reduced to F^ = 3H regardless of the dynamic 
details. Such a condition also provides a natural solution to the criticism of Yokoyama and 



Linde |]TT| for warm inflationary scenarios. 

Although an exponential warm inflation might be a viable scenario within our approach, 
probably, the more realistic one is a power law picture. It will happens, for instance, if 



the potential is proportional to the kinetic term as recently suggested in some quintessence 
cosmologies (see p^p!9| , p0| and Refs. therein). The simplest case is provided by a scalar 
field obeying the equation of state, = w p^, where w G (—1,0) is a constant parameter. 
In this case, the potential and are proportional to 0^ 

Using the new warm condition {T^ = 3H,pr = constant), the scalar field energy density is 
readily obtained from (p^. One finds 

P, = P,. [^J y [-^J , (17) 

with the quantity y defined by 



Note that if ?/ = at a finite time t* > ti, equation (|T7|) yields p^(t*) = showing that all 
the energy stored in the scalar field has been converted into photons. However, some of the 
scalar field energy must still be retained at the end of the infiationary process, in such a 
way that P(f,itf) ~ fewp^, where tf ^ t*. In fact, the end of infiation is defined by i? = 0, 
which happens for p^^. = — Pr (""^ < — |)- Naturally, the final of the infiationary stage 
is a dynamic condition, and as such, it does not means that the scalar field decaying process 
has already finished. This sort of coincidence may be only artificially adjusted so that some 
radiation entropy is presumably produced after infiation. 

Given the initial condition pr ^ p^^ (fi<^. ~ 1), as long as the scalar field dominates, the 
solution to the scale factor is 

R f3{l + w)\^^ 



Ri \ 2H, 



-1 



(t_t.)3(l+») , (19) 



showing that R{t) evolves like a power law, and a strict infiationary scenario may take place 
for ti? < — i. The corresponding number of "e- folds", = /// ^ c/t, is given by 

iV=— ^— Inf^) , (20) 

3{l + w) \tiJ ' ^ ^ 



and as shown in Table 1, it can be large enough to solve the main cosmological problems. 
In order to compute the expression for the potential V{(j)), we first consider the evolution of 



with respect to R. Using (|13|) and (|T^) one may see that the solution for the 
scalar field is 



3(l + w) , f R \ 



where Hi has been substituted by J^^P(t>i- Now, inserting ( pi]) and ( ]T9| ) into (0) we find 



(l-w) / \/247r(l + w) 



,,exp -^^ ^ ^(0-0.) , (23) 

\ '^pi I 

which shows that for w > —1, the potential V{(j)) is of the exponential type. 

Now, we recall that the ratio is usually estimated taking into account that the present 
entropy is Sp ^ Sf ^ 10®'^, which holds if after tf the universe evolves adiabatically. Thus, 
if the decaying process continues until t* > tf (p^ <ti Pr), more entropy is produced which 
means that the entropy at should be somewhat smaller than 10®"^. In what follows we 
neglect the entropy produced within this period. As one may check, if it is included, the 
upper limits derived below on the values of u are relaxed. 

By considering that the very early universe is radiation dominated, it is a good approx- 
imation to set Pr = ^gT'^ and Si = ^gT^ Ri, where g ~ 10^ is the effective number of 
massless relativistic degrees of freedom. S ince Ri — 2tj ~ 3 X 10-1° Gev'\ where U = 10"^^ 
s, T ~ 3 X 10^3 Qgy^ have pr ~ 2.7 x 10^^ Gev^ and Si ~ 3.2 x 10^^ /^From |f = (If 
one finds ^ ~ 10^^, and inserting this value into (P^, it follows that 



At ^ — ^ — - X io(36«'+2) s (24) 
3(1 + ^;) ' ^ ^ 

which yields the duration of the inflationary phase in terms of u. Recalling that when 

inflation ends the field energy density satisfy p,^^ = —j:^ Pr and expressing p,^^ as lO^pr, 

where m is an arbitrary constant, one obtains from (|lj 
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(i-'-'+a^rT^)!''-''^-'^. (25) 

which means that m > only if w < — |. 

By imposing the natural constraints p^^ < Ppi, where Ppi ~ 10^^ Gev^ is the Planck 
energy density, one obtains from the inequality (10™ + 3(1^) < 102i-72(i+«')^ which 
provides the upper limit w < —0.7 for the quintessence parameter. Therefore, one may 
conclude that the scalar field energy density in the beginning of inflation satisfies the 

required conditions, pr <^ p,/,^ < Ppi, only if w G (—1, —0.7). 

Table 1 display the limits on the main quantities in this power-law warm picture. Note 
that the prescription pr <C p^- is always satisfied. The duration of inflation is also greater 
than in the isentropic inflation. The number N of "e-folds" increases with the decreasing of 
w, and is always large enough to solve the flatness-horizon problem. 

Summarizing, in this letter we have derived a new interaction term between the decaying 
inflaton field and a thermalized bath. This coupling term is a direct consequence of ther- 
modynamic arguments, and it depends explicitly on the nature of the created particles, now 
parameterized by their enthalpy density. In this way, a consistent thermodynamic treatment 
defining clearly what we mean by a sustained thermalized bath has emerged. In particular, if 
r,^ ^ 3H, its temperature increases in the course of time regardless of the specific solution to 
R{t). We recall that such a constraint has been used as the key ingredient for warm scenar- 
ios, however, as we have demonstrated, a weakened condition (F^ ~ 3H) is not only enough 
to define warm inflation but holds regardless of the dynamic details. This means that there 
is no physical preference for exponential or power-law inflation since the evolution rate of 
the scale factor is now completely independent of the warm condition. Once the decay width 
of the scalar field and the volume expansion rate have been adjusted, the evolution of the 
Universe is exclusively defined by the energetic sector, that is, by the assumptions involving 
the relative values of V{(f)), 0^ and pr = constant. This new viewpoint has been illustrated 
by a simple warm power-law scenario, where the dynamics happens with no slow-roll regime. 
The model depends on the equation of state, = w p<^, which fixes the form of the potential 
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V{(f)) (exponential type). The duration of inflation is greater than in the warm scenario 0, 
and the model satisfies the usual constraints if the cosmic parameter lies on the interval 
— l<w< —0.7. For comparison, recent limits derived from a battery of observational tests 



(including measurements from SNe type la), constrain this parameter to he lu < —0.65 pO 



Finally, we stress that our approach can also be applied for more general models with scalar 
field coupled to a thermalized bath (generic values of F,^) . The extra bonus is the complete 
control on all the thermodynamic aspects. These thermodynamic results points to a new 
and somewhat more natural way for implementing the basic idea of warm infiation, and, in 
principle, should be justified from more fundamental calculations. 
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TABLES 



w < 


At (sec.) 


P4>i (Gev^) 


p<i, (Gev^) 


N > 


-0.71 


5.7 X 10-23 


4.9 X 1020 


1.77 Pr 


59 


-0.75 


1.2 X 10-25 


6.9 X lO^s Pr 


1.6 Pr 


61 


-0.8 


1.1 X 10-26 


2 X 10^5 


1.5 /9r 


63 


-0.9 


5.3 X 10-3° 


2.3 X 10^ 


1.2 /9r 


72 



TABLE I. Upper values for w and other basic quantities. At, p^- ( p^^) and N denote the 
duration of inflation, the energy density of the scalar field in the begin (end) of inflation, and the 
"e-fold" number, respectively. 
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